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Abstract. We use LUQAS, a sample of 27 high resolution high signal-to-noise UVES quasar 
(QSO) spectra (Kim et al. 2004), and the Croft et al. (2002) sample together with a set of 
high resolution large box size hydro-dynamical simulations run with the code GADGET-II to 
recover the linear dark matter power spectrum at z > 2 and at scales of 1 — 40 h -1 Mpc. These 
scales cannot be probed by any other observable. We address some of the uncertainties in the 
theoretical modelling of the Lyman-a forest structures such as the difference between full hydro- 
dynamical simulations and a more simplified scheme based on the HPM (Hydro-Particle-Mesh) 
technique. We combine these data sets with the Cosmic Microwave Background (CMB) data of 
WMAP in order to get tighter constraints on cosmological parameters. We focus on the recovered 
values of the power spectrum amplitude, the primordial spectral index and the running of the 
primordial spectral index. By considering models of slow-roll inflation we give constraints on 
inflationary parameters. We explore implications for the mass of a warm dark matter particle, 
including light gravitinos and sterile neutrinos, and the perspectives of constraining the neutrino 
mass. 

Keywords. Cosmology: observations, theory, large-scale structure of the universe, cosmological 
parameters - quasars: absorption lines 



1. Introduction 

The prominent absorption features blue- ward of the Lyman-a emission in the spectra 
of high-rcdshift quasars (QSOs) are now generally believed to arise from smooth density 



fluctuations of a photoionised warm intergalactic medium (e.g. Rauch (1998) I. This has 
opened up the possibility to probe the density fluctuation of matter with the flux power 
spectrum of QSO absorption lines (Croft et al. (1998); McDonald et al. (2000); Gnedin 
& Hamilton (2001); Hui et al. (2001); Croft et al. (2002); McDonald et al. (2004); Viel, 
Haehnelt & Springel (2004)). The flux power spectrum is mainly sensitive to the slope 
and amplitude of the linear dark matter power spectrum for wave-numbers in the range 
0.002 < k (s/km) < 0.05, roughly corresponding to scales 1-40 comoving Mpc/h. We 
use a suite of high-resolution hydrodynamical simulations and the flux power spectrum 
obtained from LUQAS (Large sample of UVES QSO Absorption Spectra), together with 
the published flux power spectrum of Croft et al. (2002) , to infer the linear dark matter 
power spectrum at z = 2.125 and z = 2.72. The idea is to rely on accurate high res- 
olution and large box-size hydrodynamical simulations to model the bias function b(k) 
which relates the flux to the linear dark matter power spectrum: Pp(fc) = b 2 (k) P(k). 
Moreover, we will briefly address some of the uncertainties that arise from the theoretical 
(numerical) modelling of the Lyman-a forest structures. 
We will then combine these estimates, in the framework of Monte Carlo Markov 



Chains (Lewis & Bridle (2002) |, with the Cosmic Microwave Background (CMB) results 
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Figure 1. Fractional differences between the flux power spectra computed from full hydrody- 
namical simulations and HPM simulations with different resolutions of the grid size. The shaded 
area indicates the wavenumber range used for the recovery of the underlying linear dark matter 
power spectrum. 



of WMAP ( Sp ergel et al. (2003) I to investigate the following issues: i) recovery of cosmo- 
logical parameters such as as (the power spectrum amplitude) , n s and n run (the spectral 
index and its running); ii) constraints on parameters describing slow-roll inflationary 
models; Hi) constraints on dark matter particle masses and neutrinos. 



2. The LUQAS and the Croft et al. (2002) samples 

The LUQAS sample consists of 27 QSO spectra taken with the Ultra- Violet Echelle 
Spectrograph (UVES) on VLT. Most of the spectra have been taken as part of the 
Large ESO Observing programme UVESLP (P.I.: J. Bergeron). The median redshift 
of the sample is z = 2.125 and the total redshift path is Az — 13.75. The typi- 
cal signal-to-noise ratio is ~ 50 and the pixel size is 0.05 A. For a more detailed de- 
scription of the sample and the data reduction we refer to Kim et al. (2004) We com- 
bine this sample with that of Croft et al. (2002), which consists of 30 Keck HIRES 
spectra and 23 Keck LRIS spectra and has a median redshift of z — 2.72. These 
two data sets seem to be in agreement with the flux power spectrum obtained by 
McDonald et al. (2004)[ |Viel et al. (2004)| fro m a large sample of low-resolution SDSS 
spectra. Although some differences between the different data sets are present, at the 
level of the error bars that we are going to quote, these discrepancies are not relevant. 
In all the data sets the systematic uncertainties are larger than the statistical errors 



I Kim et al. (2004)| |McDbnald et al. (2004)||viel~Haehnelt fc Springel (2004) |. 



3. Hydrodynamical simulations: HPM vs full hydro 

Presently, one of the biggest source of error is due to the not very well understood nu- 
merical modelling of the Lyman-a forest ( Viel, Haehnelt & Springel (2004) \. The Hydro- 
Particle-Mesh (HPM) method is an approximate technique to simulate the low density 
IGM ( Gnedin fc Hui (1998) I and ha s been extensively used to recover cosmologic al pa- 
rameters from Lyman-a forest data ( McDonald et al. (2004)[ Seljak et al. (2004) I. This 
method is based on the numerical implementation of an effective potential, that mimicks 
gas pressure, into a particle mesh code that simulates collisionless dark matter. We de- 
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cided to test the accuracy of this method in reproducing the flux power spectrum, which 
is the quantity that is used to infer the dark matter power spectrum, using state-of-the- 
art hydrodynamical simulations run with the GADGET-II code. The results are shown in 
figure^where we plot the fractional differences between the flux power spectrum of HPM 
simulations and the full hydrodynamical one. We can see that if the resolution (N ff id) of 
the HPM method increases then the differences can be as small as 5-10% at z — 3, 4 (mid- 
dle and left panels), in the wavenumber range of interest in the analysis (represented by 
the shaded area). At z = 2 there are strong discrepancies even at large scales that are due 
to the large amount of shock-heated gas present in the full hydrodynamical simulations 
that cannot be properly modelled with the HPM technique. The line with overplotted 
empty triangles represents the case in which the spectra have been scaled to match the 
same mean flux (both in the HPM and full hydro runs), a procedure which is usually 
adopted by simulators. We note that this scaling determines a better but still not perfect 
agreement of the two flux power spectra. Caution is thereby needed when using HPM 
simulations to infer cosmological parameters and an accurate calibration of these simu- 
lations with full hydrodynamical ones is compelling (Vicl, Hachnclt & Springcl (2005)). 



4. Results 

4.1. Cosmological parameters 
We combined the linear dark matter power spectrum as inferred by Viel, Haehnelt & 



Springel (2004) with the CMB data of WMAP ( |Spergel et al. (2003)| ). The results in 
terms of cosmological parameters are here summarized: as = 0.94 ±0.08; n s = 0.99 ±0.03 
(for models with no running, la error bars); n s = 0.959±0.036 and n run = — 0.033±0.025 
(for a model that includes running). The data is thus consistent with a scale- free power 
spectrum and no running of the spectral index. 

For slow-roll inflationary models we get the following constraints: ey = 0.032 ± 0.018, 
rjv = 0.020 ± 0.025 and £y = 0.015 ± 0.014. The constraint on the tensor to scalar ratio 



is r = 0.499 ± 0.296 (see Viel, Weller fc Haehnelt (2004) for the relation between these 



parameters and the spectral index and its running). The results above are in agreement 



with those of the SDSS collaboration (Seljak et al. (2004) I. 



4.2. Warm dark matter models 

Candidates of dark matter particles are generally classified according to their velocity dis- 
persion which defines a free-streaming length. On scales smaller than the free-streaming 
length, fluctuations in the dark matter density are erased and gravitational clustering is 
suppressed. The velocity dispersion of Cold Dark Matter (CDM) particles is by defini- 
tion so small that the corresponding free-streaming length is irrelevant for cosmological 
structure formation. That of Hot Dark Matter, e.g. light neutrinos (see following section), 
is only one or two orders of magnitude smaller than the speed of light, and smoothes 
out fluctuations in the total matter density even on galaxy cluster scales, which leads to 
strong bounds on their mass and density. Between these two limits, there exists an inter- 
mediate range of dark matter candidates generically called Warm Dark Matter (WDM). 



Viel et al. (2005) have explored two different scenarios: pure warm dark matter mod- 
els that differ from ACDM models by a suppression of power below a scale a (Figure 
2); mixed models with CDM and gravitinos (or any warm dark matter particles that 
decouples with a ~ 100, shown in Figure 3). 

We obtain a lower limit (in the case of pure warm dark matter models) of m w dm > 
550 eV (2a) for early decoupled thermal relics and m w d m > 2.0 keV (2a) for sterile 
neutrinos (that have never been in thermal equilibrium). In the case in which the gravitino 
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Figure 2. Scale of break for pure WDM thermal models. The power is suppressed below the 
scale a (Mpc/h) and the suppression depends on the mass of the wdm particle. The line shows 
ID-likelihood (upper-left panel, the marginalized result is shown as continous line) ; 2D-likelihood 
(lcr and 2a contours) in colors (the lines represent the marginalized results while the coloured 
contours arc for the mean likelihoods) for Qwdm, o"g and n 3 vs the inverse of the wdm particle 
mass. 

density is proportional to its mass we find an upper limit m grav < 16 eV (2a). This 
translates into a bound on the scale of supersymmetry breaking, A susy < 260 TeV, 
for models of supersymmetric gauge mediation in which the gravitino is the lightest 
supersymmetric particle Viel et al. (2005)| 

The values we get for the parameter a (which is related to the scale of free streaming) 
and f x = ^grav/^DM are the following: a (Mpc/h) = 0.06 ± 0.03 and f x — 0.05 ± 0.04 
(la error bars). Thereby, the concordance ACDM model is perfectly consistent with the 
data. 

4.3. Neutrino masses 

We have explored constraints on the quantity T,m L ,(eV), which is the sum of the mass of 
the light neutrino species. Both the CMB and large-scale-structure are important tracers 
of neutrino mass. At the time of decoupling neutrinos are still relativistic but become 
non relativistic later in the evolution of the universe if their mass is sufficiently high. 
As for WDM, but at larger scales, neutrinos free-stream out of their potential wells, 
erasing perturbations on smaller scales, producing a suppression of the amplitude, and 
leaving a characteristic feature at the transition scale. In Figure 01 we plot 2D likelihood 
contours for several cosmological parameters. The lines indicate the same likelihoods 
when marginalized over the full set of parameters while the coloured contours are the 
mean likelihoods. The results shown here have been obtained by using the WMAP data 



set, the 2dF galaxy power spectrum ( Percival et al. (2001)| ) and the Lyman-a constraints 
of Viel, Haehnelt & Springel (2004) We obtain the following values for the marginalized 
likelihoods: a 8 = 0.84 ± 0.08, n s = 0.98 ± 0.03, t = 0.15 ± 0.06, n m = 0.33 ± 0.06, Em„ 
(eV) =0.33 ± 0.27 (la error bars and for models which do not consider the running of 
the spectral index and tensor contribution). 
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Figure 3. Contribution of a light gravitino to the dark matter density fx = ^Igrav/^DM 
(ID-likelihood, upper left panel). The plots are the same as in figure 2 but for a model of mixed 
cold+warm dark matter. Contour plots for the cosmological parameters are now plotted vs the 
gravitino mass (or any particles that decouples with a g, — 100) . 



5. Conclusions 

The Lyman-a forest, the absorption by neutral hydrogen seen in the spectra of back- 
ground QSOs, is a powerful laboratory for cosmology. We have used two samples of 
high resolution high signal-to-noise QSO spectra to infer the linear dark matter power 
spectrum at z > 2 and at scales not probed by any other observables. The present un- 
certainties in the inferred values of the cosmological parameters are mainly due to the 
systematic errors involved in the numerical modelling of the Lyman-a forest data with 
hydrodynamical simulations. We have explored the differences between full hydrodynam- 
ical simulation and an approximate scheme based on the HPM approximation. 

We have then combined the small scales estimates of the Lyman-a forest with larger 
scales data sets such as the CMB and the 2dF galaxy power spectrum. The error bars on 
cosmological parameters that can be inferred from the combined analysis of CMB and 
Lyman-a forest data improve by a factor two those obtained by the CMB data alone. 
The final results point to a value of a§ = 0.94 ±0.08 and a spectral index n s = 0.99 ±0.03 
(lcr error bars) and no evidence for a running of the spectral index. These results are 



consistent with those obtained by the SDSS collaboration Seljak et al. (2004) We have 
set constraints on the mass of a warm (thermal) dark matter particle to be m w dm > 550 
eV and on the light gravitino to be m grav < 16 eV (2ct error bars). Further constraints 
on the sum of the three neutrinos can be obtained and we get: Smj, (eV) =0.33 ± 0.27 
(la error bar) by combining WMAP, 2dF and Lyman-a . 

Acknowledgements 

I thank my collaborators: Julien Lesgourgues, Martin Haehnelt, Sabino Matarrese, 
Antonio Riotto, Volker Springel and Jochen Weller. I thank IAU for the award of a 
grant. The simulations were run on the COSMOS supercomputer at the Department of 
Applied Mathematics and Theoretical Physics in Cambridge. 



6 



Matteo Viel 




References 

Croft R. A. C, Weinberg D. H., Katz N., Hernquist L., 1998, ApJ, 495, 44 

Croft R.A.C., Weinberg D. H., Bolte M., Buries S., Hernquist L., Katz N.,Kirkman D., Tytler 

D., 2002, ApJ, 581, 20 
Gnedin N.Y., Hamilton A.J.S., 2002, MNRAS, 334, 107 
Gnedin N.Y., Hui 1998, MNRAS, 296, 44 

Hui L., Buries S., Seljak U., Rutledge R. E., Magnier E., Tytler D., 2001, ApJ, 552, 15 

Kim T.-S., Viel M., Haehnelt M.G., Carswell R.F., Cristiani S., 2004, MNRAS, 347, 355 (K04) 

A.Lewis and S.Bridie, 2002, Phys. Rev. D 66, 103511 

McDonald P., Miralda-Escude J., Rauch M., Sargent W.L., Barlow T.A., Cen R., Ostriker J.P., 

2000, ApJ, 543, 1 

McDonald P. et al., 2004, submitted to ApJ, |astro-pn /0405013 

Percival et al. 2001, MNRAS, 327, 1297 

Rauch M., 1998, ARA &A, 36, 267 

Seljak U., et al., 2004, |astro-ph/0407372] 

Spergel D. N. et al., 2003, ApJS, 148, 175 

Springel V., Yoshida N., White S.D.M., 2001, NewA, 6, 79 

Viel M., Haehnelt M.G., Springel V., 2004, MNRAS, 354, 684 

Viel M., Weller J., Haehnelt M.G., 2004, MNRAS, 355, 23L 

Viel M., Haehnelt M.G., Springel V., 2005, astro-ph/0504641 

Viel et al., 2004, proceedings of the XXXIXth Rencontres de Moriond on "Exploring the Uni- 
verse", La Thuile, (Italy), astro-ph/0405584 
Viel M., Lesgourgues J., Haehnelt M.G., Matarrese S., Riotto A., 2005, PRD, 71, 063534 



